Quantized vortices in ^He droplets: a quantum Monte Carlo study 
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We present a diffusion Monte Carlo study of a vortex line excitation attached to the center of a 
''He droplet at zero temperature. The vortex energy is estimated for droplets of increasing number 
of atoms, from TV = 70 up to 300 showing a monotonous increase with A'^. The evolution of the core 
radius and its associated energy, the core energy, is also studied as a function of A'^. The core radius 
is ~ 1 A in the center and increases when approaching the droplet surface; the core energy per unit 
volume stabilizes at a value 2.8 Ka~^ {a = 2.556 A) for > 200. 

PACS numbers: 67.55.-s,36.40.-c,02.70.Ss 



Quantized vortices are one of the most outstanding sig- 
natures of superfluidity. They have been widely observed 
in liquid '*He at ternperatures below the critical temper- 
ature Tx = 2.17 K More recently, vortices have also 
been detected in dilute Bose [3| and Fermi [3j gases when 
they are magnetically trapped, their observation being 
considered the most clear indication of the achievement 
of their superfluid phases. '^He droplets produced in free 
jet expansion experiments [3| are expected to be also su- 
perfluid due to their very low temperature (T = 0.38 K). 
However, in this finite system the search for a direct sig- 
nature of their superfluid character is more difficult. On 
the one hand, indirect evidence of superfluidity has been 
driven from the determination of the rotation spectrum 
of molecules adsorbed into them [5] . This very interesting 
phenomenon has been considered the microscopic version 
of the famous Andronikashvili experiment. On the other 
hand, the detection of vortices in droplets would be an 
even more conclusive proof of their superfluidity. Never- 
theless, no signal of straight or circular vortices has been 
yet observed in experiments with *He droplets, in spite 
of some theoretical arguments in favor of their possible 
metastability. 

The stability of vortex excitations in ''He droplets was 
first studied by Bauer et al. 6]. They concluded that 
vortices are not stable in droplets due to the high exci- 
tation energy required, compared to the usual temper- 
ature at which they are produced. Later on, Lehmann 
and Schmied studying cold droplets, smaller than 
the ones analyzed in Ref. [1] were led to a different 
conclusion. They stated that in small droplets, where 
only surface excitations (ripplons) are relevant, vortices 
should be stable against decay. Finally, density func- 
tional (DF) calculations ^,'9'] have predicted that below a 
critical atomic number a linear vortex pinned to a dopant 
atom or molecule can become stable, with a lifetime long 
enough to allow for its experimental detection. 

The excitation energy of a vortex line in a ^He droplet 
as a function of its number of atoms N is one of its more 
fundamental properties. An accurate calculation of its 
value is crucial to elucidate its possible stability and for- 



mation probability in jet expansions. To our knowledge, 
there is only one previous microscopic calculation of the 
excitation energy associated to a vortex in a droplet. 
This calculation, performed with the path integral Monte 
Carlo (PIMC) method [13], was carried out for a iV = 500 
^He droplet and the energy obtained was more than a fac- 
tor two smaller than DF predictions In the present 
work, we present zero-temperature diffusion Monte Carlo 
(DMC) results of vortex energy and vortex structure in 
■^He droplets for different number of atoms. 

A vortex is an excited state of the TV-particle Hamil- 
tonian which corresponds to an eigenstate of the angular 
momentum operator. Actually, it is an eigenstate of the 
z component of the angular momentum, L^, with eigen- 
value hNl where I — 1,2,... is the quantum of circula- 
tion. In a general form, the imaginary-time dependent 
wave function of a vortex in a quantum liquid can be 
written as 



(1) 



By imposing that '^^ is an eigenstate of the angular mo- 
mentum Lz, Feynman ll| obtained his famous proposal 
r!(R,i) = ;(/)f(R), with (/.f(R) = Ej=i,iv^j' being 
the j-th polar coordinate (in cylindric coordinates). 

With the decomposition ([1]), the Schrodinger equation 
for 5'u(R, <) splits in two coupled equations, one for the 
modulus 

- ? = D{VVtf'^ - Z3V2$ + (F(R) - E)* , (2) 
ot 

and one for the phase 
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with D — h /2m. If a fixed form for the phase Q is as- 
sumed, the equation for the modulus, Eq. becomes 
the usual Schrodinger equation with one additional term 
(the first one on the right hand side). In this approxi- 
mation, known as fixed phase (FP), the vortex acts like 
a static external potential. Using Feynman's expression 
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for the phase, = IJ^iLi this potential results 



N 



(4) 



with Pi the polar coordinate of particle i in cylindric co- 
ordinates. Therefore, in the FP approximation the prob- 
lem of having a vortex inside the droplet is reduced to a 
different Hamiltonian {H = H + Vv) in the Schrodinger 
equation to be stochastically solved. This approach was 
used in the PIMC calculation of the vortex energy in a 
''He droplet and also recently in a Monte Carlo study 
of the excitation energy of vortices in trapped diluted 
Bose gases [12 1. 

In the present work, we will assume that the vortex 
line is fixed in the z direction of the center-of-mass (CM) 
reference system of the droplet. The vortex state is then 
an eigenstate of the L^cm operator, accounting for the 
translational invariance of the Hamiltonian. In this case, 
the resulting potential (V^) that must be added to the 
Hamiltonian in the droplet geometry is 
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Coordinates in Eq. ([5|) and hereafter are referred to the 
center of mass. The potential Vv is very similar to the one 
for the bulk (H]), but now written using CM coordinates, 
and a small correction of order 1 /N is introduced. 

The FP approach to the excitation energy of a vor- 
tex line with the Feynman's phase could be thought 
as a too crude approximation for an accurate micro- 
scopic treatment of the problem. An a priori better 
method would be to consider the superposition of clock- 
wise and anticlockwise vortices, which are degenerate in 
energy, and using the resulting wave function (orthog- 
onal to the ground state) as a guiding wave function 
in DMC. This leads to a Fermi- like problem due to the 
non-positivity of the excited wave function which can be 
approximately solved in the fixed-node (FN) approxima- 
tion. This method was used in Ref. [l3| for studying a 
vortex in a two-dimensional geometry. However, the re- 
sults there obtained showed that FN and FP predictions 
are almost compatible. Moreover, both Ref. [l3| and Ref. 
14j studied the possible improvements upon Feynman's 
phase by introducing backflow correlations on it and con- 
cluded that their effect on the excitation energy is very 
small (< 1% in the energy per particle). Therefore, it 
is sound to consider that the FP method is also a good 
enough approximation to describe the vortex attached to 
a droplet. 

We have carried out DMC simulations of ^He droplets 
with number of atoms TV = 70, 128, 200, and 300. To 
extract the excitation energy associated to the vortex line 
we have calculated the energy of the droplets for both 
the ground state and the droplet with a fixed vortex in 



the center of mass. The ground-state trial wave function 
used for importance sampling in the DMC method is a 
Jastrow form. 
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The first term in the square bracket of Eq. ([5]) is a McMil- 
lan correlation factor accounting for dynamical correla- 
tions induced by the interatomic potential V{rij), and 
the second one is a gaussian correlation to take into ac- 
count the finite size of the droplet. 

When the vortex is present, the Hamiltonian changes 
due to the potential induced by the FP approximation 
([5]). Therefore, we introduce an additional correlation 
factor f{p) to reduce the variance of the energy estima- 
tion, 
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$.(R,i) = $o(R)[]/(p,) , 
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with pi the radial position of particle i. Among the differ- 
ent options for /(p), discussed for the bulk in Ref. [l3|, 
we have chosen 



f{p) ^ 1 - e-p''^ 
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with a a variational parameter related to the radius of 
the vortex core. We have taken a = 1 A, value which 
coincides with the old estimation of the core radius by 
Rayfield and Reif [l5j and the more recent one by Sadd 
et al. [l3|. This function satisfies /(p) when p 
0, reflecting the repulsive character of the potential Vy, 
and approaches 1 far from the core. We have checked 
that the explicit form of the function /(p) satisfying both 
boundary conditions is not really important and that, as 
expected, the energy of the system does not depend on 
it. 

The other two variational parameters b and a appear- 
ing in the trial wave function ^ have been optimized 
by means of variational Monte Carlo (VMC) calcula- 
tions. The optimal values are b = 1.19(t and a = 0.6a~^ 
(cr = 2.556 A) and their dependence with the number of 
atoms N is negligible in the range studied. 

The main purpose of this work, i.e., the determination 
of the excitation energy associated to a vortex line at- 
tached to the CM of the ^He droplet, has been studied 
by other groups using mainly density functional theory 
0, H, 0] . Dalfovo et al. [1] proposed a liquid-drop formula 
for the energy dependence on the number of particles. 



(9) 



A, P and 7 being parameters which were obtained by 
fitting Eq. ^ to their DF results. The final set of pa- 
rameters was A = 2.868 K, /3 = 1.445 K and 7 = 0.313 
K. 
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N 


70 


128 


200 


300 


Eo[K] 
E,[K] 
E4K] 


-223.6(2) 
-188.0(2) 
35.7(3) 


-497.4(5) 
-446.6(4) 
50.8(6) 


-861(1) 
-799(1) 
61(2) 


-1376(2) 
-1302(1) 
74(2) 



TABLE I: Total energies of "^He droplets with (Ei) and with- 
out (Eq) a vortex inside, and vortex excitation energy AU 
energies are in units of Kelvin; in parenthesis, the statistical 
errors. 



The specific dependence on TV in Eq. ^ is derived 
using the hollow-core model. In this model, the local 
vortex energy is integrated over all the volume occupied 
by particles 

E^N)^ [ dV^ , (10) 
Jv P 

resulting in 
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Notice that one can rewrite Eq. pTjl in terms of N just 

using R = roN^/'^, where = (3/(47rZ?o))^^'^i arriving in 
such a way to the same dependence on N as Eq. 
However, the hollow-core method is a too simple approx- 
imation to identify the parameters A, f3 and 7 using, for 
example, a constant density Dq. 

Calculations performed using DF Q and MC evalu- 
ate the vortex energy as the difference between the total 
energies of the droplet with and without (ground-state) 
a vortex {Ey{N) = Ei{N) - Eo{N)). In Table I, we 
present our DMC results of the energies Ei{N), Eq{N), 
and E^{N) for '^He droplets with N = 70, 128,200, and 
300. The vortex energy E^ comes from the difference of 
two energies which increase with the number of atoms N 
and therefore its statistical error also increases with N. 

The results contained in Table I show that the exci- 
tation energy Ey increases monotonously with N. This 
can be more clearly observed in Fig. [T] where the present 
DMC results are compared with the liquid drop formula 
dl]) as reported by Dalfovo et al. [1]. One can observe 
that the DF estimation reproduces better our results for 
the larger droplets than for the smaller ones. This is 
what one a priori expects since DF approximations work 
better for large droplets where application of mean field 
theory is more justified; in spite of this, the difference is 
only ^ 4% for the smallest droplet studied [N — 70). 

The repulsive potential induced by the vortex creates a 
hole when p — > with a characteristic modulation of the 
density. This behavior is shown in Fig. 2 where we have 
plotted the density profiles D{p) for different values of z 
along the vortex axis. The slice at z = corresponds to 
the one in the center of the droplet and is the density pro- 
file with the highest peak. When z increases, the radius 
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FIG. 1: Vortex energies for ''He droplets with A'^ = 
70, 128, 200 and 300. The energy curve as a function of N 
corresponds to the DF results from Ref. @]. 
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FIG. 2: Density profiles for a A = 200 *He droplet with 
a vortex inside. They are represented as a function of the 
radial coordinate p and for different slices corresponding to 
different z values on the vortex axis. From top to bottom 
the curves correspond to increasing values of z from z = in 
the center of the droplet up to z values close to the surface. 
For comparison, we also show the radial density profile of the 
droplet without the vortex (solid line). 

of the slice is progressively smaller and also the oscilla- 
tions are depressed as z approaches the droplet radius. 
In the center of the vortex, the density goes to zero in 
agreement with the DMC results derived en Ref. [l3| for 
homogeneous 2D liquid ''He. A relevant parameter in the 
microscopic description of a vortex is the size of its core, 
what is called the core radius ^. There is not a single def- 
inition of ^ but different definitions lead to quite similar 
results. We have used the criterion of considering for ^ 
the position of the maximum in the azimuthal circulating 
current Jq (p) . If the vortex is described by means of the 
Feynman approximation, as in the present work, it fol- 
lows that Je{p) — D{p)/ p. Therefore, we can get a direct 
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FIG. 3: Local energy Ey[p) ior a N = 70 ""He droplet with a 
vortex inside. From top to bottom the curves correspond to 
increasing values of z from z = in the center of the droplet 
up to z values close to the surface. 



estimation of ^ from the density profiles D{p) shown in 
Fig. 2. In the center of the droplet f = 1.0 ± 0.1 A and 
it increases vifhen z approaches the value of the droplet 
radius up to ^ = 1.4 ± 0.1 A. 

In a homogeneous system, and at large distances from 
the vortex axis, the excitation vortex energy is usually 
decomposed in the form 



(12) 



with A a constant, Ec the core energy, and p the radial 
coordinate in cylindrical coordinates. This model, which 
reproduces very well DMC results of the vortex energy in 
a 2D geometry 13[ , splits the vortex energy in a constant 
term associated to the hole around the vortex axis and 
a logarithmic term containing the hydrodynamic tail. In 
the case of droplets, the behavior of the excitation energy 
is more complex due to their inhomogeneity along the z 
vortex axis and their finite size. In the function Ey{p) 
one finds, superimposed to the monotonously increasing 
law ([T2|) . a decaying trend to zero when the surface of 
the droplet is reached. This behavior is shown in Fig. 3 
where we have plotted the function Ey{p) for a = 70 
''He droplet. The function Ey{p) shows clearly a peak 
corresponding to the core of the vortex, especially for the 
inner slice, and a decay to zero in the surface. The size 
of the droplets studied is too small to see the signature 
of the hydrodynamic tail but we can estimate the energy 
of the core Ec by summing the local energy up to the 
estimated vortex core radius ^. Extending this sum along 
all the vortex axis we obtain an energy which increases 
with N: Ec = 9.74(4) K and 19.2(3) K for = 70 
and 300, respectively. Normalizing Ec with respect to 
the volume of the core, the core energy per volume unit 



approaches a constant value of 2.8 Ka ^ for droplets with 
N > 200. 

Summarizing, we have carried out the first microscopic 
calculation of the properties of a vortex line attached to a 
^He droplet at zero temperature. The energies obtained 
are in good agreement with DF estimations and therefore 
give additional confidence on their predictions. More- 
over, the magnitude of the core radius and core energy 
has also been studied for the first time in these inhomoge- 
neous system. The core radius is ~ 1 A in the center and 
increases when approaching the droplet surface; the core 
energy per unit volume stabilizes at a value 2.8 Ka~^ 
for A^ > 200. We hope our results can help and stimu- 
late further experimental work towards the long-standing 
purpose of detecting vortex excitations in ^He droplets. 
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